Supercoil formation in DNA denaturation 
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We generalize the Poland-Scheraga (PS) model to the case of a circular DNA, taking into account 
the twisting of the two strains around each other. Guided by recent single-molecule experiments on 
DNA strands, we assume that the torsional stress induced by denaturation enforces formation of 
supercoils whose writhe absorbs the linking number expelled by the loops. Our model predicts that, 
when the entropy parameter of a loop satisfies c < 2, denaturation transition does not take place. 
On the other hand for c > 2 a first-order denaturation transition is consistent with our model and 
may take place in the actual system, as in the case with no supercoils. These results are in contrast 
with other treatments of circular DNA melting where denaturation is assumed to be accompanied 
by an increase in twist rather than writhe on the bound segments. 

PACS numbers: 87.15.Zg, 36.20.Ey 



Thermal denaturation of double stranded DNA 1] has 
been of recent interest due to its relevance to protein syn- 
thesis, polymerase chain reaction (PCR) and microarray 
technologies. The denaturation transition has been ex- 
tensively studied theoretically mainly by means of two 
models: (a) the Poland-Scheraga (PS) model Q which 
considers the opposite bases to be either bound with a 
certain energy gain or unbound (when part of a ssDNA 
loop) and ignores the twisting of the strands around each 
other; and (b) the Peyrard-Bishop model [3|] which is a 
one dimensional model in which complementary bases in- 
teract by a distance dependent potential. Both models 
have been used to study, e.g., the nature of the melt- 
ing transition & & 0, @| and the bubble dynamics of 
the DNA [1, 0,11,13 • 

It has recently been shown that 
the PS model predicts a first-order melting transition if 
one properly takes into account the self-avoidance of the 
chains [lEEEl- 

The original PS model treats the DNA as a long lad- 
der, without considering the twisting of the two strains 
around each other. It can be argued that this feature 
is irrelevant for the thermodynamics of an open ended 
chain, since the twisting strain can be released by the 
rotation of the chain ends. However, this assumption is 
no longer appropriate if the chain ends are not free to 
rotate or alternatively, for circular DNAs such as plas- 
mids in bacteria. In this case, upon heating up to the 
melting temperature the two strands can no more fully 
depart from each other, since the chemical bonds that 
assemble the sugar backbone are still intact. We assume 
topoisomcrases and other topology modifying agents are 
not present in the solution. 

Then, the partition function is restricted to a sum over 



configurations with a fixed linking number (the number 
of times one strand rotates around the other). A well- 
known theorem 12] states that the linking number (LN) 
is the sum of twist and writhe, where twist refers to the 
sum of the subsequent stacking angles along the DNA 
and writhe is associated with the geometry of the DNA's 
center line and measures the amount of twist absorbed 
by the excursions of the backbone. 

Past attempts to include DNA's helicity in the PS 
model have considered a twisting strain associated with 
the modified stacking angle of subsequent base pairs ac- 
cumulating upon loop formation. This assumption leads 
to the conclusion that the transition either changes its 
nature (becomes of higher order) or disappears alto- 
gether, depending on the treatment of the self-avoidance 
effects |l3l | and the applied external torsion [3] ■ In both 
works, as well as here, it is assumed that the LN accom- 
modated by the loops is insignificant. The reason is that 
the single strands may be considered as random chains 
and their winding angle is rather small [l5| . 

Recent experiments on single DNA chains, however, 
point to a different possible mechanism for the absorp- 
tion of LN expelled from the denaturated loops. Exper- 
imental measurements of the torsional response of long 
strands under fixed stretching force F s show that when 
F s is small, the chain almost immediately undergoes a 
buckling transition, forming a supercoil that absorbs the 
externally introduced LN [16( . This mechanism of har- 
boring the LN in the modified conformation of the center 
line (writhe) rather than in augmented basepair stacking 
angles is a familiar phenomenon frequently observed in 
telephone chords. The presence of a similar phenomenon 
in DNA calls for a re-examination of thermal denatu- 
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ration in a setting where the LN associated with the 
loops is transferred to supercoils formed by locally re- 
laxed DNA segments. Indeed, recent experiments point 
to the presence of such a mechanism during plasmid dc- 
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FIG. 1: A depiction of a microscopic configuration of the 
generalized Poland-Scheraga model discussed in the text. The 
superscripts (d),(c) and (s) refer to denaturation loop, coil and 
supercoil regions, respectively. 

In Fig. [1] we depict the generalized PS model we con- 
sider for this purpose. In the model, a particular base 
pair may be either unbound (in a loop), bound in a coil, 
or bound in a supercoil. Unlike in the original PS model, 
a looped configuration is allowed only if the rest of the 
chain can be rearranged to form the required amount of 
supercoils for the conservation of the total LN. Let the 
linking number stored in the natural twist of a relaxed 
DNA in I subsequent base-pairs be accommodated in the 
writhe of a supercoil segment of total length £' = 81. We 
will set 5 — 1 below, since a general treatment, although 
possible, does not bring new insight. Furthermore, it is 
reasonable to expect that loops form on the coils but not 
on the supercoils which are relatively rigid structures. 
As a result, formation or expansion of a loop on a coil is 
accompanied by an equal increase in the total length of 
supercoils in the system. 

Let Eb < be the binding energy of bound pairs in the 
coil or supercoil state and E s > be the cost of increasing 
the length of a supercoil by one bp (e.g., due to the bend- 
ing rigidity). For simplicity, we assume the minimum 
size of a supercoil to be one base pair and one expects 
E s < —Ef, for the DNA. The corresponding Boltzmann 
factors are u> = exp[— j3Eb] and v = exp[— (3{Eb + E s )], 
where /3 = 1/kT. The Boltzmann factor corresponding 
to the configuration in Fig. [1] is then, 
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Here if represents the length, in units of base pairs, of 
the loop i. The segment of bound pairs separating the i 
and i + 1 loops is composed of alternating subsegments 
of coiled regions of length Z?- and supercoiled regions of 
length Ifj . The entropic contribution of a loop of length 
Z is f2(2Z) = As l /l c , where c is the universal entropic 



parameter of a loop, s is a non-universal constant, and 
A is a parameter which incorporates the cooperativity 
parameter, i.e. the Boltzmann weight associates with 
the initiation of a loop. Typically A ~ 10 -4 and it is 
weakly temperature dependent [19(. Note that, some 
inter-loop regions may accommodate several supercoils, 
whereas some may have none. 

Let Ld, Lc and L$ be the total length of the denat- 
urated regions, the coiled and the supercoiled regions, 
respectively, in a given configuration. The length of the 
molecule is L = Ld + Lc + Lg. The canonical partition 
function is a sum over the contributions of all microscopic 
configurations with total DNA length L, and with a fixed 
LN, namely with Ld — L$. In the grand-canonical en- 
semble, the two constraints are relaxed by introducing 
two fugacities /x and z, which contribute an extra weight 
z l^l d -l s |- each configuration. The two fugacities are 
determined by taking the appropriate derivative of the 
grand partition sum Q: 
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Ignoring end effects which contribute terms of order L 
to the partition sum, one finds that the grand partition 
function Q can be expressed as 



Q{z,n) = V{z,p) + V(z,ii)U(zii)V(z,ii) 
= V/(1-UV) , 
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where, following Ref. 

U(z[i) = 

with, V(z) 



En(2/)(z/x)' = A$ c (sz/i) , (4) 
i=i 

V(z)/[l-W(z/ l i)V(z)}, (5) 

oo 

=u)z/(l-uz) , 
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w( X ) = j2(^y = Y 



(6) 
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The functions U , V , and W represent the grand sums for 
a loop, a coil, and a supercoil, respectively, and <& c (x) is 
the Polylog function. Note that, the functional form of 
the grand sum in Eq. ([3]) is similar to that of the original 
Poland-Scheraga model, except that the "propagator" 
for the coil regions is now dressed to accommodate 
an arbitrary number of supercoils. The price paid for 
conserving the linking number is that its associated 
fugacity [i needs to be calculated as a function of z at 
each temperature. 
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After some algebra Eq. Q reduces to the more transpar- 
ent relation dw{ d z J^ + = 0. Using Eqs.© and 
©, we obtain the following transcendental equation for 
fi(z): 
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In order to study the nature of the denaturation tran- 
sition (when it exists) we consider the thermodynamic 
limit (L — > oo) by focusing on the relevant pole of Q(z) 
at U(z*n*) = l/V(z*/n*), where ^* = n{z*) through 
Eq.©. Substituting l/V = l/V - W, we obtain 
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= A$ c (sz*n(z*)) . (8) 



The average bound pair density in coiled and supercoiled 
segments is given by 
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and 6e 



dlogz* 
dlogv 



(9) 



respectively. Therefore, a phase transition is associated 
with a singularity in the temperature dependence of the 
solution z* of Eqs.© and ©. 

Before proceeding, let us consider the simpler picture 
where the DNA chain is not circular, but, nevertheless, 
can form supercoils. Assuming the ends are bound and 
free to rotate, we set /i = 1 in the grand sum, since the 
system is now insensitive to the linking number. Eq.© 
alone suffices to describe the phase transition in this case. 
Let the LHS and the RHS of Eq.© with (j, = 1, be named 
F(z) and G(z), respectively. Then, F(z) is a smooth, 
monotonically decreasing function of z for < z < 1/v, 
with F(0 + ) — +oo. Similarly, G(z) is a smooth, mono- 
tonically increasing function of z for < z < l/s (and 
divergent for z > l/s), with G(0) = 0. A phase tran- 
sition exists if the smallest z* that satisfies Eq.© ex- 
hibits a singularity as a function of temperature. This 
is the case only if c > 1 (so that & c (sz) remains finite 
as sz — > 1). In addition one requires F(l/s) > G{l/s) 
at infinite temperature where u) = v — 1. Equivalently, 
after substitution. 



s-l/(s-l) > l + AQ , 



(10) 



where £ c = $ c (l)- Given a suitable set of the phenomeno- 
logical constants that satisfy these conditions, the model 
exhibits a phase transition of first (second) order for c > 2 
(1 < c < 2). This mechanism is described in detail in 

In. 

Having established the existence of a phase transition 
in the unrestricted case and with the possibility of su- 
percoils, let us now turn to the circular DNA, where /x is 
determined by Eq. © in order to ensure LN conservation. 
For a circular DNA, two regimes emerge: 

For 1 < c < 2, a solution of Eqs.© and © with 
szfi = 1 does not exist, since Eq.© with $ c _i(l) = oo 



dictates fx(z) = vz, whereas in Eq.©, $ c (l) < oo. In 
fact, 1/V(z,[i) — > — oo at szfi = 1 for all temperatures, 
ensuring a smooth variation of function of tem- 

perature (Fig. [2]). Hence, the second-order melting tran- 
sition found for the unrestricted DNA is absent when the 
linking number is conserved. 
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FIG. 2: Absence of a melting transition for c < 2. The pole 
of the grand sum varies smoothly at all temperatures, since 
l/V —* — oo at szfi = 1. The shown plot of l/V is qualita- 
tively same for all temperatures. 

For c > 2, substituting sz*fi* = 1 in Eq.© and 
Eq.©, and picking the smallest positive solution for z* , 
one obtains 



C_ and 
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with C± = y/1 + l/(4AC c -i) ± a/1/(4AC c -i). The criti- 
cal temperature is given by 
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Since v 1 ! 2 /u> > 1, a phase transition exists only if s is 
sufficiently large. In particular, if 



s > ACc-i 



(13) 



when A < 1. The transition temperature is reduced 
relative to the unrestricted case roughly by a factor 
- log s (s/AC c -i}. With s w e 12 5 as used by the MELT- 
SIM scheme [20( one expects a transition to take place 
at a finite temperature (although the values of s and A 
optimized for our model may be different). The two qual- 
itatively different regimes separated by the boundary in 
([15]) are depicted in Fig.©. 

We conclude that the circular DNA undergoes a denat- 
uration transition only if c > 2. Performing a variational 
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analysis near sz*[i* — 1, it is straightforward to show 
that the phase transition is of first order for all values 
of c > 2. Also note that, the fraction of bound pairs in 
coiled and supercoiled segments, 8> and Of, in the high 
temperature phase follow from Eqs.© and (jTTJ) as 

9> = and 6> = 1/2 . (14) 

The high temperature phase, although devoid of coiled 
regions, is not fully unbound. Half of the chain is in a 
supercoiled state while the remaining half unbinds, max- 
imizing the amount of entropically favored loops. 
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FIG. 3: Two melting scenarios in the model with c = 2.115, 
s — 4.68, E s — 0, and (a) A = 0.1 with a first-order transition 
at u>crit — 4.18 as found from Eq. (|12p . (b) A — 1.0 with no 
phase transition. Full circles are the simultaneous solutions of 
Eqs, (17l8[ ) corresponding to the thermodynamic limit. These 
parameter values which are relevant for self-avoiding walks 
on a cubic lattice [llj suitably demonstrate the two possible 
scenarios, but are not intended to fit the experiments directly. 

These results contrast with the solution of the PS 
model without twist Q , as well as with the earlier calcu- 
lations on the effect of locked-in twist on DNA melting 
transition [l3, 14\. The no- twist PS model predicts a sec- 



ond order denaturation transition for 1 < c < 2, whereas 
if supercoil formation is taken into account, this second- 
order transition is absent when the LN is conserved. Ear- 
lier attempts to incorporate twist, where the loop forma- 
tion is penalized by the overtwisting of the bound seg- 
ments, found that the first-order transition which exists 
for c > 2 in the original PS model becomes of higher or- 
der. In the proposed supercoiling scenario, the transition 
for c > 2 remains first order. 

Our findings should also apply to long DNA chains, 
where the twist expulsion through the ends may be hin- 
dered by kinetic effects (see, e.g.. [ill]). A more general 
framework where both supercoiling and twisting effects 
are incorporated is called for, in order to confirm the free 
energetic preference for the proposed mechanism over the 
alternative overtwisting scenarios. 
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